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ABSTRACT 

In the context of massive fragmenting protoplanetary discs, iBossI (J1998D suggested 
that grains can grow and sediment inside giant planet embryos formed at i? ~ 5 AU 
away from the star. Several authors since then criticised the suggestion. Convection 
may prevent grain sedimentation, and the embryos cannot even form so close to the 
parent star as cooling is too inefBcient at the s e dist ances. Here we reconsider the 
grain sedimentation process suggested by iBosJ ()1998l ) but inside an embryo formed, 
as expected in the light of the cooling constraints, at i? ~ 100 AU. Such embryos are 
much less dense and are also cooler. We make analytical estimates of the process and 
also perform simple spherically symmetric radiation hydrodynamics simulations to test 
these ideas. We find that convection in our models does not become important before 
a somewhat massive (~ an Earth mass, this is clarified in a followup paper) solid core 
is built. Turbulent mixing slows down dust sedimentation but is overwhelmed by grain 
sedimentation when the latter grow to a centimetres size. The minimum time required 
for dust sedimentation to occur is a few thousand years, and is a strong function of 
the embryo's mass, dust content and opacity. An approximate analytical criterion is 
given to delineate conditio ns in which a gian t embryo contracts and heats up faster 
than dust can sediment. As iBoss et al.l ()2002[ ). we argue that core formation through 
grain sedimentation inside the giant planet embryos may yield an unexplored route 
to form giant gas and giant ice planets. The present model also stands at the basis of 
paper III, where we study the possibility of forming terrestrial planet cores by tidal 
disruption and photoevaporation of the planetary envelope. 
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1 INTRODUCTION 



The Gravitatio nal Instability (GI) model for giant planet 
formation (e.g.. lBodenheimeilll974l : lBosslll997h is one of the 
two w ell known models for planet formation (e.g.. lWetheril]| 
[l990|). The model has been thought to be unable to pro- 
duce solid cores observed in giant plane ts in t he So l ar Sys - 
tem (e.g.. JFortnev fc Nettelrnannl 120091 ') until IBossI (| 19981 ): 
IBoss et al.l ( 2003 ) demonstrated that dust growth and sed- 
imentation may realistically occ u r ins ide gaseous proto- 
planets. However, IWuchterl et al.l (|200(J ) argued that giant 
planet embryos may become convective within the first ~ 
100 years of their ev olution, and that co n vection may inhibit 
:rain sedimentation. [Helled et all (|2008l ): [Helled fc SchubertI 
2008 ) have recently confirmed that convection may indeed 
be a serious obstacle to solid core formation inside the giant 
planet embryos. T hese authors obtained much smaller cores 
than the model of IBossI (|l99g ) predicts. More recently, the 
GI model for giant planet formation has been criticised on 
the grounds that the cooling time at the inner, e.g., i? ^ 10 
AU disc, is not short enough to permit a continuous grav- 



itational contrac tion of gaseous clumps (e.g.. lRafikovll2005l : 
iRice et al.|[2005[ ). 

Whilst this criticism of JBossJ ([1993) ideas appears to 
be relevant, there are now extrasolar giant planets observed 
at distances out to hundreds of AU f rom their parent star , 
where the core accretion model (e.g.. IWeidenschilling|ll98(]| : 
IWetheriUlllQQOJ ) for planet formation may not have been able 
to create solid cores before the gaseous disc dissipates. It 
seems rather likely that the GI model is needed to explain 
these planets. Simulations of large (e.g., tens of AU) discs 
that employ realistic cooling in the optically thick regime do 
show formation of gaseous clump s with masses of a few to a 
few tens of Jupiter mass es (e.g., IStamatellos fc WhitworthI 
I2OO8I : iMeru fc Bate|[2()lCll ). which presumably cool and con- 
tract into giant planets or brown dwarfs. 

It is timely to reconsider the gravitational disc insta- 
bility model taking into account the updated constraints on 
the disc fragmentation (e.g., iRafikov 2005|). To investigate 
this complicated and non-linear problem properly, we divide 
presentation of our work into several papers - (1, this paper) 
before and (2, Nayakshin 2010b, to be submitted) after the 
solid core formation, and (3, Nayakshin 2010c, submitted to 
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MNRAS letters) interplay of the processes occurring inside 
the embryo with external influences of the parent disc and 
the star. 

The goal of this paper (referred in later papers as "paper 
I" ) is to study dust growth and sedimentation inside the gi- 
ant embryo but excluding the solid core formation. We limit 
our study to isolated non-rotating giant planet embryos or 
first cores, the terms that we use inter-exchangeably below. 
While the study of such a limited problem may appear too 
academical, we shall use these results extensively in papers 
II and III. There is also a possibility that dust sedimenta- 
tion inside first cores may also be relevant to the process of 
low mass star formation. Therefore we hope that analytical 
results presented here may be useful for future researchers 
in more than one field. 

The gist of the "new" disc fragmentation constraints is 
that the giant planet embryo mus t star t off m uch less dense 
than previously assumed. Namely, iBosg (| 19981 ) assumed that 
the protoplanet has a mean density of order '^ 10~* g cm"''. 
Such high a density is needed to overcome the tidal fields 
from the parent star at several AU distances. However, at 



100 AU, the initial embryo density of order 



10" 



cm is expected (see paper II for detail). Optical dep ths of 
such embryos are much less t h an of those s t udied by iBosd 
||1998D : IWuchterlet"ai] (|2000D : iHelled et all (|2008D . There- 
fore we find that the giant planet embryo remains radiative 
long enough to permit dust growth and sedimentation. In 
fact, convection becomes important only after the solid core 
assembly as the energy release by the solid core rises to as 
much as 10^* - lO^^'erg s"^ (paper II). 

The paper is organised as following. In iJ21 we intro- 
duce the so-called "first cores" and argue that the initial 
state of giant planet embryos should be similar to that of 
the first cores. In Sj3]we study the internal evolution of iso- 
lated first cores, defining the maximum gas accretion rate 
onto the first core under which the core can be considered 
isolated. We then attempt to carry out the calculations ana- 
lytically as far as we can, building a simple analytical model 
to follow the cooling and contraction of the first cores, and 
grain growth and sedimentation inside them. We also con- 
sider the influence of turbulent mixing and grain shattering 
via high speed collisions that may delay the onset of grain 
sedimentation. 

In later parts of the paper, ^ and ^ we use a simple 
spherically symmetric radiation hydrodynamics code with 
a two-fluid approach to capture grain growth and dynam- 
ics to allow for a simultaneous evolution of the gas and the 
dust. We flnd higher opacity, larger metalicity and lower 
mass first cores to be the most promising sites of grain sedi- 
mentation. Turbulent mixing mainly slows down rather than 
forbids grain sedimentation; although in corners of the pa- 
rameter space it may be crucially important. We conclude 
that whether the grains sediment inside the first core or not 
mainly depends on whether it contracts and hence heats to 
the grain vaporisation temperature faster than the grains 
can grow. Finally, in the Discussion section we overview the 
main results of our work, compare those to previous work, 
and consider implications for the field of planet formation. 
We argue that the planet formation community should take 
a detailed look at the giant instability disc model supple- 
mented by a realistic internal evolution for the giant embryo, 
and also embryo dynamics in the parent disc. 



2 FIRST CORES AND GIANT PLANET 
EMBRYOS 

2.1 First cores in studies of star formation 

Stars form from larger gas reservoirs due to gravitational 
contraction of the latter. The simplest case of a spherical, 
uniform, non-rotating, a nd initially at rest, constant den- 
sity cloud is well known (|Larsonlll969l '). The collapsing gas 
is initially isothermal, as cooling is very effective. However, 
as the innermost region of the collapsing gas fragment be- 
comes denser, the free fall ti me there becomes sh orter than 
the cooling time of the gas (JMasunaga fc Inutsuka J_99y) . 
The region switches from an isothermal to an adiabatic be- 
haviour. Heated up, the gas is able to stop the further col- 
lapse by thermal pressure forces. The central part of the 
region settles into a hydrostatic balance. This region is log- 
ically named "the first core". 

The minim um mass of the first cores is given b y the 
"opacity limit" (|Reeslll976l : ILow fc Lvnden-Belllll976l) of a 
few Jupiter masses (although iMasunaga fc Inutsukal 1 19991 , 
note that the term is somewhat physically misleading), and 
the size is ~ 10^** cm. The conditions in a typical non rotat- 
ing Solar mass molecular gas cloud puts the first cores in a 
situation where the mass supply rate is continuous and rela- 
tively high (~ a few xlO~^ M© yr~^). In this case new mat- 
ter lands on the first core so fast that it does not have time to 
cool, and therefore the behaviour of the gas nearly adiabatic 
(lLarsonlll969l : iMasunaga et al.lll998l : IMasunaga fc Inutsukal 
|2000|) . They contract and heat up only due to more and more 
mass piling up on them from the parent gas cloud. When 
their mass exceeds about 50 Jupiter masses, the cores reach 
the temperature of about 20007^ and the gas density about 
10~* g cm~^. Dissociation of molecular hydrogen then en- 
sues. This provides an additional energy sink and initiates 
a "second collapse" which finally turns the first core into a 
true proto-star as compact as a few Solar radii. 

To understand why the first cores do not collapse fur- 
ther until they reach a mass of about Mmax = .05 Mq 
JMasunaga et al.l Il998l : IMasunaga fc Inutsukal bOOd l , con- 
sider the virial temperature of the first core. As will be shown 
latter (equation [T^]), an adiabatically evolving first core has 
the virial temperature of about 
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150 K 



Affc 



0.01 M, 



4/3 



(1) 



where Affc and Rf^ are the mass and the radius of the first 
core, respectively. The core becomes hot enough for hydro- 
gen to disassociate only when it reaches the mass of about 
50 Jupiter masses. 



2.2 Giant planet embryos are similar to first cores 

Consider now the properties of giant planet embryos at the 
moment of their formation in a protoplanetary disc. In order 
to fragment gravitationally, a massive gas disc must satisfy 
two criteria. First of all, i t must be massive/dense enough. 
In particular, the iToomrel ([196J) Q-parameter must be less 
than unity: 
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where Cs is the gas isothermal sound speed; Q, = 
[GMt,/ RpY''^ is Kepler's angular frequency; M, is the mass 
of the star, assumed to be larger than that of the disc, and 
E = 2Hp is the surface density of the disc. Here p is the 
disc vertically averaged density, and H = CgSl^^ is the disc 
vertical scale-height. Rearranging equation[21 we see that at 
(5 = 1 the disc mean density is equal to the tidal density of 
the star. 
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(3) 



to within a factor of order unity, where R2 = Rp/IQO AU. 

The second condition for disc fragmentation is the con- 
dition on the rate of cooling. Expressed i n terms of the 
cooling time tco oi, the requirement states (|Gammiel I2OOII : 
iRice et al.|[2005l ): 

tcooi ;S 3!^-' . (4) 

This implies that the critical radiative cooling rate per unit 
gram of disc material is 



Ac 
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f'cool 



ciQ. 



3(7-1) 



(5) 



where 7 is the specific heats ratio for gas. At the same time, 
th e compressional heating rate for gravitational contraction 
is (|Masunaga fc Inutsukall 19991 ) 



Agrav = cl {'iTlGpf'^ = 2^''^ clO, 



(6) 



where we have used equation[3]to eliminate p. It is now clear 
that Acrit ~ Agrav at the point of disc fragmentatiorij. This 
is intuitively correct, as otherwise the disc could heat up 
faster than it could fragment, increasing Cs and returning 
itself back to stability (e.g., Q > 1). 

We have shown above that a fragmenting self- 
gravitating gas disc satisfies Acrit ~ Agrav. This is also the 
condition that sets the " opaci t y limit for fragmentation" 
in star formation (iReej Il976l: ILow fc Lvndcn-Belll 1 19761 : 
iMasunaga et al.l 119981 : iMasunaga fc In utsuka 1999)), which 
yields a minimum mass for fragmentation of a few Jupiter 
masses. This is of course not a coincidence as in both cases 
gravitational collapse is stopped if the gas is unable to con- 
tinue cooling rapidly enough. 

Therefore the properties of the gaseous clumps in the 
disc at their inception should be similar to that of the first 
cores of a few Jupiter masses. As with the first cores, the 
clumps may gain more mass from the disc. However, the 
giant planet embryos in the disc have to compete for gas 
with other embryos, as the disc is likely to fragment on 
many rather than one clump. Collisions between clumps 
may unbind the clumps partially or completely. Further- 
more, it is well kn own that above the "transition" mass 
Mt -^ 2M:^{H/Rf iBate et all 120031 ') the planet opens up 
a gap in the disc due to gravitational torques between the 
planet and the disc, and that the gap opening c urtails gas ac- 
cretion onto th e planet by orders of magnitude (jLubow et al.l 
ll999l : lBate et al.ii2003 ). Giant planet embryos with mass of 



^ This statement is accurate within a factor of order unity, e.g., 
the same order of accura cy to which the vertically integrated 
accretion disc equations jShakura &: SunvaevI 119731 ) are to be 
trusted. 



the order of 10 Jupiter masses are in this gap-opening regime 
as long as H/R ^ 0.2 for a Solar mass star. 

Therefore, for simplicity, we study the evolution of 
planet embryos at a constant mass in this paper. The mass 
of the gaseous clump is a free parameter below. We also 
use the terms "first core" and "giant planet embryo" inter- 
exchangeably throughout. 



3 INTERNAL EVOLUTION OF FIRST CORES 

3.1 Physical properties of adiabatic first cores 

IMasunaga et al.l (| 19981 ) show that the structure of the first 
cores can be well approximated by a polytropic sphere, 
with the polytropic constant fixed by the gas temperature 
and the density pad at which the gas switches from the 
isothermal to the adiabatic behaviour. The typical values 
of these are often taken in the literature as Tinit = 10 K 
and Pad = 10"^"^ g cm"''. However, pad does depend on the 
opaci ty of the material and the ambient gas temperature 
Tinit IMasunaga fc Inutsukalll999h . Following these authors, 
we take opacity to be dominated by dust, in the functional 
form 



n{T) 



Ko 



©"■ 



(7) 



where kq is the opacity at T = 10 K, and reasonable values 
of a are thought to be between 1 and 2. With this, the 
adiabatic density threshold is best described by the following 
(JMasunaga fc Inutsukall 19991 ): 



Pad = 5 X 10 
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T, 



g cm 



(8) 



whe re k« = kq/O.OI and Ti = Tinit/IOK. 

IMasunaga et al.l (|l998|) give expression for Rtc as a func- 
tion of the central density. Eliminating the density in favour 
of the sphere's mass and radius, we obtain 



-Rfc 



17.5 AU m^^^^TipZli^ 



(9) 



where p_i3 = lO^^Pad, Ti = Tlnit/lO and mi = 
Affc/0.01 M0. Eliminating the critical density by using equa- 
tion (8] we get 



Ric ^ 6.0 AU m^'-''^T; 



4/9 



(10) 



Note that the temperature dependence in this expression is 
at most linear, and the opacity enters in an even weaker 
power. The initial size of the first core is hence varies rather 
little. The mean density of the first core is defined by 

l+4o 



K, 



4/3 



(11) 



SMfc in-12 -3 2rrn 

Pmean = „;; = 6.6 X 10 g Cm niiT^ 

and is a very strong function of the mass of the first core. 
The temperature of the first core is of the order of the virial 
temperature, which we estimate to be 



_ 1 GMi^p 

J- vir — -, 



146 K mj/^Tj" 



l+4o 
9 



-4/9 



(12) 



3 fcs-Rf, 

From t his it f ollo ws that the sound speed in the first core is 

c, = ^fcT/p ^ y/GMfJR{,, 



Cs = 0.7 km s m-^ T^ 



-2/9 



(13) 



The column depth of the first core is 
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^ Mic 3 _2 5/3rri^^Tr^ -8/9 

Efc = — ^2" = 1.18 X 10 g cm jTij^ ^1 '^j 
and the optical depth of the first core is 



(14) 



K(rviOSfc = 11.8 X (14.6)"mi ^^ MTi,k*) , (15) 



where (f}i {T\ , k, ) = T-^ 



K,f 



From this expression 



it is obvious that the first cores are always optically thick 
The radiative luminosity of the first core is 

AnRlaBTX- 



AcLfc 

or numerically, 



(16) 



(17) 



1+4q/9 



Lfc = 2.2 X 10''(14.6)""mi ^ (^2(Ti,«:*) , 

where function (f>2{Ti,Kf) — j'^ ~ -"^ °' >' f^~ 
cumbersome looking power-law function which does not vary 
much for a between 1 and 2, however. For example, for a = 
1, 02(21, K,) = Tj^' k7 , whereas for a = 2, (t>2{Ti,K,) = 
Tj K^ . Nevertheless, it does become quite inconvenient 
to write down expressions for a general value of a, and so 
we shall present the results only for either a = 1 or a = 2 
below. For convenience of future reference, 



1.5 X W^'m'/'T'/'k:'^" 



for a = 1 
for a — 2 



(18) 



Note that the luminosity is lower for a = 2 than for the 
other case. This is due to a higher opacity of the first core, 
e.g., the larger exponent in the opacity law. The binding 
energy of the polytropic sphere, 

3 GMl 



Ebin 



1.8 X lO"*^ erg m^^T^ 



-4/9 



5 Ru 
allows us to calculate the cooling time of the sphere as 

-Bbind 



Lcc 



(19) 



(20) 



For the smaller value, a = 1, 
icooi = 380 years i-nl^'''T~'^^'''K,l^^ 
whereas for a — 2, 
tcooi = 5700 years mtT^ 



(21) 



(22) 



Another useful time scale is the dynamical time of the 
core, which is of the order of the free fall time for the core: 



idyn = (Gpmcan) ^ ~ 50 yrs m^ ^T^ ** hil 



2/3 



(23) 



Note that the cooling time is indeed longer than dynamical 
time, as should be for a hydrostatic gas configuration. 

3.2 Gas-starved regime for the first cores 

Having reviewed the adiabatic evolution of the first cores, 
we can now estimate when cooling of the cores becomes im- 
portant. As the core gains (or perhaps looses mass), the 
adiabatic temperature time derivative is 



rfTvi 



dt ) ad 9A/fc dt 3 dt 



(24) 



where we used eauation ll2l On the other hand, the temper- 
ature derivative due to radiative cooling is 



dTs,\ 
dt 



J- vir 

i'cool 



(25) 



(see also equation[21]below). We define a critical "radiative" 
first core accretion or decretion rate as the one for which 
|(drvir/di)rad| = {dTvii / dt) a_d ■ FoT luass gaiu or loss rate 
below this critical rate, the radiative cooling is the dominant 
driver of the cloud's thermal evolution. From equation 1251 
and 1241 we find 



dAffc 



dt 



rad 



3Mfc 
4t 



cool 



2 X lO"-" 
1.3 X 10" 



Moyr 
-^Msyr 



1 1/3 

m 



-.4/3,-1/9 



19/9^1/3 



-Irril 



for a = I 
for a — 2 



(26) 



(27) 



Note further that as the first cores contract, their cooling 
time tcooi increases ( i|3.3l below). and hence the critical radia- 
tive accretion rate actually drops with time from the value 
found in equation 1271 

We can now compare the critical radiative accretion 
rates with the accretion rates of the first cores in the con- 
text of a Solar mass gas cloud collapsing, for which the first 
cores accumu l ate mass at the rate dMjc/dt ~ 10~^ Mq yr~^ 
(Larsonlll969l:lMasunaga et al.lll998l: iMasunaga fc Inutsukal 



|2000|). These rates are higher than the critical radiative ac- 
cretion rates, and therefore the first cores should indeed 
evolve quasi-adiabatically in these conditions. 

In the opposite case, |dMfc/df | ^ (dAffc/df)rad, the evo- 
lution of the first cores will be governed mainly by their ra- 
diative cooling. This is the case that we henceforth study in 
this paper. 

3.3 Cooling of isolated cores 

Staying within the analytical approach, we calculate the 
cooling evolution of the first cores, assuming that they re- 
main polytropic spheres. The solutions obtained in this way 
will be later compared to a numerical calculation that does 
not use this assumption. 

We obtain the solution by solving the equation 
d-Ebind/di = I/fc, where the plus sign is used because we 
defined the binding energy to be positive. For convenience, 
we introduce dimensionless temperature T(t) of the first core 
as the ratio of the current virial temperature T{t) to the ini- 
tial virial temperature of the core at formation, as given by 
equation 1121 We also introduce a constant for a given Affc, 
to = icooi(O), which is the cooling time of the first core at 
formation, given by equation 1201 

Equation [TS] shows that luminosity of the first core 
changes with time as 

Lfc (X RiT^-°- (X f -" ; (28) 

where we used the fact that Rt^T oc is a constant for a given 
Mfc. Since the binding energy scales as -Ebind « T{t) at a 
fixed Affc , the equation for the evolution of the dimensionless 
temperature of the first core can be written as 



dT _T 



to 



which is solved as 

m = 

Using this, we find for the cooling time evolution 



l + (l + a)- 
to. 



(29) 



(30) 
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tcool(t)=tcool(0) + (l + a)i 



(31) 



which in the limit t ^ fcooi(O) shows that the cooling time is 
always 1 + a times the current time. Note that this implies 
that the first core's cooling becomes more inefficient as it 
contracts. The luminosity of the core evolves according to 



Ltc(i) = ifc(0) 



l + (l + a) 



(32) 



where Z/fc(0) is the initial luminosity of the core, given by 
equation 1 181 Note that equation[32]is an increasing function 
of time, so that the first cores do become brighter as they 
contract, but the binding energy of the core increases with 
t even faster, and this is why the cooling time is a growing 
function of time as well. 

We can now define the time scale ivap defined as the 
time needed for the first core to heat up to the temperature 



Tvap at which grains vaporise (see i|4.2.2l T^-, 
1400-?^, depending on the size of the grains) 
is obtained by solving 



-^ van — -^ vir 



1 



to 



1200 - 
This time scale 



(33) 



where Tvir is the initial virial temperature of the first core as 
given by equation 1121 As an example, in the limiting cases 
when tvap ^ to, the results are: 
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(35) 



for a = 2. These time scales are very long, and are strongly 
dependent on the mass of the first core. The term neglected 
in the last two equations becomes important for more mas- 
sive cores that are already hot at their birth. They may reach 
the vaporisation temperature very quickly. 

Figure [T] shows the evolution of the virial temperature 
(simply referred to as temperature of the first cores here- 
after) for several values of Mfc as a function of time for the 
two limiting opacity cases. The curves are computed using 
equation I3UI The upper panel presents the faster contract- 
ing case a = \ and the lower one corresponds to a = 2. In 
the case of the lower opacity (q = 1), the relatively massive 
cores, Mfc ^ 20A'fj, contract quite rapidly during the first 
few thousand years. In terms of grain growth and sedimen- 
tation, these cores are not very promising, as we shall see 
below. In contrast, for M^ ^ lOMj for q = 1, and for all 
values of Mfc for a = 2, the cores take t > 10* — 10^ years 
to contract to the grain vaporisation temperature. Grain 
growth is plausible in these cases. 

We also plot first core's mean densities as a function of 
time in Figure [2] As in Figure [l] the curves are terminated 
when the core's temperature rises above 1400ii', as then the 
grains would be vaporised. We see that the lower mass first 
cores evolve (contract) significantly before vaporisation of 
grains occur. This again hints at less massive first cores as 
the more promising sites of grain growth and sedimentation, 
as they simply stay cooler for longer. 



Figure 1. Virial temperature as a function of time in the ana- 
lytieal model for cores of different masses, A/fc = 3, 5, 10, 20 and 
30 Mj, from bottom to top curves, respectively. The top panel 
is for opacity power-law index a = 1, whereas the bottom panel 
is for a = 2. The plots are terminated at T = 1400 K, when even 
the largest grains would vaporise rapidly. Note that lower mass 
cores could support grain growth for longer as they are initially 
cooler. 



3.4 



Time scales for grain growth and 
sedimentation 



IBossI (|l998l ) studied formation of gaseous giant planets by 
gravitational instability in a proto-planetary disc. He con- 
sidered a gaseous clump that would e ventually con tract to 
form a proto-planet of mass M — IMj . IBossI (| 19981 ) demon- 
strated that if the conventional arguments for dust growth 
and sedimentation are specialised for the spherical geometry 
of the clump, then one can expect a rather rapid dust growth 
and sedimentation which may culminate in a formation of a 
heavy elements core. 

Because of the specifics of the gravitational disc insta- 
bility model for planet formation, the gas density of 10~^ 
g cm^'' was used for the gaseous clump. As shown in ij3.ll 
our gaseous clumps are much less dense and are also cooler 
in th e ir init ial stages. Nevertheless, the line of arguments of 
IBossI (|l998l ) for the dust growth and sedimentation model 
can be simply rescaled to the problem at hand. Our treat- 
ment thus follows his model. 

We shall first assume that turbulent motions in the first 



6 S. Nayakshin 











y y 




10« 


r 




y y 
yy 












F 








y^ y" y 


i 


10^ 


r 




. -/.^x^ - 


>, 




















U) 








-^ .- y^ '' / 


0) 


in-^° 


^ 






a 




? ~ 




^y^ y -' - 


CO 






CD 

2 


10" 


r- 




-''' / : 
















^--- 




^■ 




10^^ 




-^■ 


1 



10^ 



lo-" 

time [years] 



10" 



10" 




10^ 



10^ 10" 

time [years] 



10' 



Figure 2. Same as in Figurc[T] but now showing the mean density 
of the first cores. The mean density curves are terminated when 
first core's temperature reaches T = 1400 K. Note that maximum 



mean densities are a few xlO 
cores but can reach ~ 10~* g cm" 



: cm ^ for the more massive 
for the least massive cores. 



core can be neglected. This allows us to make the most op- 
timistic estimates of grain growth and sedimentation. The 
importance of turbulence will be elaborated in i|3.7l and fur- 
ther tested with numerical models in §3. 

Further, we note that realistically, grains of different 
sizes are present in the cloud at any given time, due to frag- 
mentation of larger grains by high speed collisions. Only 
a fraction of grains is large enough to start sedimenting 
down. Therefore the mass fraction fg in our model should be 
thought to represent only those larger grains. A sizable frac- 
tion of the original grains is assumed to remain small and 
tightly bound to the gas. Accordingly, a constant gas opac- 
ity of the form given by equation [7| is used despite allowing 

larger grains t o sediment. 

FoUowing lSafronovl (| 19691 ) ; IWeidenschilJiii^ (|l980l ) ; iBosd 

1 19981 ) ■ we assume that t he grains can grow by the hit-and- 
stick mechanism (but see lBlum fc Wur m 20081). For a con- 
stant density core, the gravitational acceleration, Ogr, 



47r 

"gr = ——-GpfcR , 



(36) 



where R ^ Rf^ is the radial position of the grain inside the 
first core. 



There is no pressure gradient force for the grains, but 
there is a gas-grain drag force if the grain velocity, Ua, is 
different from that of gas, u. Note that Ua is always subsonic 
with respect to the gas speed of sound, as the largest velocity 
that the grain can attain is the free-fall velocity which is of 
the same order as Ca by definition. 

The drag force on a spherical body depends on the 
Reynolds number, Re — a\Au\/\Cs, where Au = Ua — u, 
and A is the mean free path for hydrogen molecules in the 
gas. The latter is relatively large, e.g., A = l/(n(7H2) ~ 40 
cm pZiQ, where p_io = 10^'^pfc is the dimensionless density 
of the first core. As we shall see below, grains will usually 
satisfy a ^ A , in which case the Epstein drag law applies, 
and one has (|Bosslll998h 



/dua_\ 
V dt J, 



picCs 



(Ua 



(37) 



drag 

Here a, Ea = paa and pa ^ 1 g cm"'' are the radius, the 
column and the mass density of the grain, respectively. 

For a larger body, a > A, the Stokes drag law applies. 
The law also depends on the relative velocity Au through the 
Reynolds number (see I Weidenschillinej|l977l ) . In the limit of 
the small Re < 1, The Stokes law yields 



f dua \ 



ipicCs A 



(Ua 



(38) 



drag ^^a 

Neglecting the complicated behaviour of the drag coeffi- 
cient for intermediate values of Re, we combine both regimes 
described by equations [37] and [38] into one, approximately, 
as 

'dUa\ _ picCs A 

IT ^ ~ 



fdUa\ 



-{Ua — u) 



(39) 



' drag ^a (1 -\~ \ 

Given our simple one size dust model, a better gas drag 
treatment appears to be excessive, but future more detailed 
calculations should utilise more careful drag force treat- 
ments. 

The equation of motion for the grain is then 

dUa PfcCs A , , 47r . 

~dr 



-{Ua 



-GpicR 



(40) 



Ea a + A 3 

Small grains, a <^ picR/pa, quickly reach their terminal ve- 
locity, so that dua/dt ~ 0, and the grains slip through the 
gas in the direction of the centre of the cloud with the sedi- 
mentation velocity 

4TTGT.aR A -fa 



"sod 



{Ua 



3cs 



A 



(41) 



Note that this velocity is proportional to R. For this reason, 
for a given radius of the grain a, grains starting at different 
R will fall to the centre at the same time. We refer to this 
time scale as the sedimentation or dust settling time: 



ts, 



R 



3cs 



A 



47rGEa a + X 



5x10^ yrs mi^(paa) 



.(42) 



where the density of the grain pa and its size a are in cgs 
units, and we assumed a <g A limit in the last step. 

For microscopic grains, say 1 micron, the sedimentation 
time scale is prohibitively long. Therefore grains must be- 
come larger before any sedimentation takes place. The early 
growth of microscopic grains is dominated by Brownian mo - 
tions of the smallest grains (|Dullemond fc Dominikll2005l ). 
Since the sedimentation velocity is a function of the grain 
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size a, grains of different size move with differential speeds. 
Hence the larger grains sweep smaller grains, leading to the 



growth of the grain's mass, rua 
drua 



(47r/3)paa , at the rate 



dt 



■na fgpUscd 



(43) 



where fg = 0.01 is the mass fraction of grains, so that pg — 
fgP is the density of grain material. This translates into the 
growth rate for a of 

da fgp 



dt Apa 



"'scd 



(44) 



One finds that this differential settling grain growth is faster 
than that due to Brownian motion for larger grains. We now 
define the grain size e-folding time scale, fc ~ a/{da/dt), 

3cs 



fe(0) 



(45) 



TVfgPfcGR 

Using R ~ Rfc, the grain growth time scale is defined as the 
time needed to increase the grain size from an initial value 
ao to the final size a: 

3cs 



U 



nfgpfcGRfc 



In- 



ao 



(46) 



Choosing a = 10 cm and ao = 10 * cm as an example. 



ln(a/ao) ~ 23. Noting that Rfc/cs — tdy 



[Gp) 



-1/2 



see that for fg = 0.01, some 2 x 10 dynamical times need 
to pass before the grains attain a size at which they can 
sediment rapidly. Using the relations from i]3.1l we have 



V(0) = 5.4 X lO^mrVra'fc'/'rr''"'"'"^^''^^^^!^ , (4T) 

in years, where /_2 = /g/O.Ol. In this equation we empha- 
sised the fact that the grain growth time estimate is obtained 
for t = 0. 



3.5 Time-dependent grain growth 

Comparing the grain growth time scale with the cooling 
times of the first cores (equations [21] and I22p . and with 
the longer vaporisation times (equations [34] and [35]), one 
is tempted to think that grain growth is too slow to be of 
importance, especially for the lower opacity index, a = 1. 
However one needs to be more careful here as the first core 
contracts with time. Therefore grain growth actually accel- 
erates as the gas cloud cools, and especially rapidly for the 
case Of = 1. 

Using equation[46]and our model for the evolution of the 
first cores, we plot the time-dependent estimate for the grain 
growth time scale in Figure [3] for the same first core masses 
as in Figures [T] and [2] Although the lower mass cases (upper 
curves in Figure [3]) initially have longer grain growth time 
scales, they are also initially cooler and stay sufficiently cool 
to avoid the second collapse for longer. Therefore, these cores 
should be promising sites of grain growth after a few hundred 
years, when they cross the line shown with diamonds, which 
is simply the current time. 

We can improve the analytical estimate for the time 
needed for grain growth (equation I47p by utilising the re- 
sults from i]3.3l We observe that igr oc Cs/{p{cRfc) oc 
(Mfc/iifc)^'^^(i?fc/Mfc) oc f(t)-^''^ at a constant M^. There- 
fore, the equation for the grain growth in a contracting first 
core is 
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Figure 3. Grain growtli time scale, tgr (equation I46I I. for time- 
evolving first cores of different masses and the two opacity cases 
(cf. Figures [T] and [2J . The diamonds show the line igr = i. Any- 
where where the igr curve drops below the diamond line, the 
grain growth is possible. The curves are terminated when the 
cores become hot enough to vaporise the grains. Lower opacity 
(a = 1) and higher mass first cores cannot support substantial 
grain growth and sedimentation. 



da 



to(0) 



1 + (1 + q)- 
lo. 



1.5 

l + n 



(48) 



This equation can be integrated analytically, and inverted 
to give a more accurate estimate for the grain growth time 
from size ao to a: 



ig 



to 



1 + a 



1 + ^^(0)^^ 



1 



(49) 



where ^ = (l-|-a)/(2.5 + a). In practically all the interesting 
parameter space, igr(O) 3> io, and the above equation can 
be further simplified to 



ig 



(2.5 + a)« 1 



igr(0)« 



(50) 



H-a 

This results in considerably shorter grain growth times. In 
particular, for a = 1, omitting the dependence on Ti, 

tg.(a) = 3.3 X 1^'mf^'r-rkV'^ i^^^^Y^ ' (^D 
and for a = 2, 
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tgr(a) = 2.3 X lO^'/Ia' fc 



1/3 fMa/ao) 



20 



2/3 



(52) 



Both of these cases show that grain growth is possible within 
a reasonable range of parameters, as the vaporisation time 
fvap can be comfortably longer than this. Furthermore, nu- 
merical experiments (Q show that grain growth time is 
about a factor of two shorter than the analytical model pre- 
dicts. 



3.6 Sedimentation for a constant grain size 

Based on our numerical experiments below, grains do not 
sediment strongly until their size reaches at least a few cm. 
Approximately then, we divide the process of grain growth 
and sedimentation into the growth phase, studied above, and 
the sedimentation phase. Assuming that grain growth slows 
down in the latter phase (due, e.g., to a too high relative ve- 
locity between the grains, see i ]3.8p . we consider the grains 
to have a constant size whilst they sediment. This simplified 
model allows analytical calculations to be carried through. 
We shall continue to approximate the first core with the 
constant density p = pfc for R ^ Rfc and zero density out- 
side, as in tj3.1l Furthermore, we shall fix these, essentially 
assuming that sedimentation occurs on a time scale shorter 
than the current cooling time of the first core. 



3.6.1 Homologous contraction in the gas-dominated phase 

We start off with the grain density radial distribution fol- 
lowing the gas distribution, thus pg{R,0) = fgptc, and zero 
grain velocity. The gas is assumed stationary. One can show 
that the grains achieve the terminal velocity u — —Us^d 
feg nation 141 1) very quickly when \va\ <^ Cg. Here we kept 
the dependence of the drag force on the grain size in the 
more general form (equation [39]) . Note that the sedimenta- 
tion time scale tscd is independent of R in the limit a ^ \ 
(see eauation l42p . Radial position of grains evolves according 
to 



dR _ R 

dt tscd 

which is trivially solved, 
t 



R{t) = Ro exp 



fs, 



(53) 



(54) 



where Ro is the initial value of 7?. R{t) is the Lagrangian 
coordinate for a grain mass shell Mg{Ro). This equation 
shows that the grains contract homologously, so that the 
grain density is independent of radius within -Rg(i) = 
-Rfc exp(— t/iscd), and is zero outside. The grain density pro- 
file keeps its top-hat shape but becomes more compact. The 
grain density increases with time as 

r 3t 1 
Ps(i) = /sPfcexp - — . (55) 



3.6.2 Grain- dominated core phase 

The equation 21] for the grain settling velocity neglects grav- 
ity from grains, which is appropriate in the initial stages of 
the process. When the density of the grains approaches and 



then exceeds that of the gas, the equation obviously be- 
comes inaccurate. The contracting cloud of grains becomes 
self-gravitating at the time 



iscif = is 



ln(/,-') 



— i.OO tged 



The grain sphere's outermost radius at that time is 



Rsclt = Rtcfg 



0.2 



0.01 



Ric 



(56) 



(57) 



Thus to follow the contraction of the grain sphere after time 
tacit, we modify the equation for settling velocity by writing 

dR _ GMg{Ro)T.a{a + X) 
dt 



R^ 



PfcACs 



(58) 



which now completely neglects the mass of the gas interior 
to radius R, as pg ^ pf^, asymptotically. Since Mg{R) = 
Mg{Ro), i.e., constant in Lagrangian coordinates, the above 
equation is solved as 



R{t) = R{ta. 



t-ts 



U 



1/3 



(59) 



which is a homologous contraction again, albeit at a different 
- accelerated by the self-gravity - rate. The density evolution 
follows the form 



Pa{i) = Ptc 



t~t. 



is, 



(60) 



where we utilised the fact that the grain density at the ini- 
tial time when the solution 1601 becomes applicable, t — tscd 
is equal to pfc. Within this simple model, all the grains col- 
lect to the centre of the first core, reaching formally infinite 
densities, at time 



is, 



Coo — isclf ~r ■ 



for fg = 0.01. 



He/fo) 



tscd — 1.87 isod 



(61) 



3.6.3 Bound "grain cluster" phase 

Before contraction of the grains into a point occurs, another 
important milestone is reached when the grain sphere be- 
comes not only self-gravitating but also gravitationally self- 
bound. In the self-contracting phase studied above, the grav- 
itational force acting on the grains is dominated by the grain 
density. Hence the grains are self-gravitating in that sense. 
However, if the outer gaseous envelope were removed in the 
beginning of that phase, gas inside the grain sphere would 
create a significant pressure gradient that could unbind the 
contracting grain-gas mix. In contrast, at later time, when 
the grain sphere contracts even further, its density is high 
enough that removal of the gaseous envelope would not un- 
bind the grains. We refer to this stage as the "grain cluster" 
one, in analogy to a star cluster. 

Consider then the question of how compact the grain 
sphere of mass fgMfc should be in order to be self-bound, 
i.e., so that the gas pressure could not unbind the sphere. 
We follow the usual order of magnitude arguments with 
which the Jean's mass can be derived. The gas pressure 
gradient is ~ pknT / p.R, whereas the gravitational accel- 
eration is {GMenc{R) / R^){p + Pg). The total enclosed mass 
is Menc{R) = (An /3){p -\- Pg)R^ . This defines the radial scale 
(Jean's length and also the size of the grain cluster) 



Grain sedimentation in giant embryos 9 



Rp^ 



3 ksT p 
A-K fiG pI 



1/2 



(62) 



where we explicitly assumed p ^ pg. The Jeans mass is 
equal to Mj = Mgc ~ {4tv / 3) pgR'^^, which can be re-written 
as 



AL. 



1/2 

AnpJ y pG 



o \ I Kb! \ IP 

Pa 



(63) 



The first two factors on the right hand side of the equation 
1631 is the mass of the first core, Mf^. Further, Mgc = fgAdfc, 
and hence we can solve for the grain density in the grain 
cluster: 



-(—1/2 

PgC ~ Pfc Jg 

Using this result in equation |j 
ter is found to be 



(64) 
the radius of the grain clus- 



Rp- 



Rfc /, 



1/2 



O.lRt, 



0.01 



1/2 



(65) 



This equation shows that the grain sphere becomes self- 
bound when it contracts to size 0.17?fc, e.g., 0.1 — 1 AU 
for a realistic range on parameters. 

Note that the grain density in the grain cluster is typi- 
cally a factor of ten or more higher than the density of the 
first core itself, i.e., as high as lO"'' — 10~^ g cm"''. 

Presumably, once the "grain cluster" phase is reached, 
nothing keeps the grain-dominated region from a self- 
gravitational collapse, so that a solid core could be formed 
l|Bosa 1 19981 1. We delay the study of this issue till a future 
paper. 



3.7 Turbulent mixing 

Turbulent mixing (|Fromang fc Papaloizoull2006l ) is a process 
in which grains are dragged along turbulent motions of the 
gas. As initially grains represent a small fraction of the total 
mass, they are essentially a trace population, and therefore 
gas turbulence drives diffusion of grains. This process hence 
tends to erase grain density inhomogeneities, opposing grav- 
itational settling of the grainS; 

Numerical simulations (jFromang fc Papaloizoul (|2006r ). 
and references thereafter) show that a simple diffusion equa- 
tion approach describes the effects of turbulence on dust 
well, provided that the diffusion coefficient, D, is chosen 
right. The diffusion equation for grains in the spherical ge- 
ometry can be written as 



dpg 
dt 



D d 
R^dR 



r2 d{Pg/p) 



dR 



(66) 



where p and pg are the gas and the grain (dust) density, 
respectively. Note that we assumed that D is a constant in- 
side the giant embryo, e.g., independent of R. The diffusion 
coefficient in case of disc can be parameterised in the form 
D = OdCa H, reminiscent of the stan dard disc viscosity pre- 
scription (jShakura fc Sunvaevlll973l ). where a^ < 1 is the 
viscosity coefficient, and H is the disc vertical scale height. 
We shall use the same approach, except that the scale-height 
H should be replaced by the radius of the first core, _Rfc: 



D = adCsRfc 



(67) 



The diffusion coefficient ad is unknown here. In the disc 
geometry, differential rotation velocity develops between 
the dust layer and the gas, driving instabilities (see, e.g., 
iGaraud fc Lin 2004'). The s imula tions of the turbulent disc 
bv lFromang fc Papaloizorj (|2006l ') showed that Od ^ 0.004, 



but there is no clear reason why their results should trans- 
late to the case we study. Therefore we shall treat ad as a 
free parameter of the model and consider the implications 
of turbulent diffusion in the "small" and "large" ad cases. 

To continue with our analytical modelling for now, we 
shall again assume the "top hat" density profile for the gas, 
and hence set p to a constant. Also, we shall combine the 
diffusion equation for dust (grains) with the mass continuity 
equation; 



dpg 
dt 



D d 
R^dR 



R- 



■ dpg 
OR 



wm [""'■•''l 



(68) 



where V — — Uscd is the grain terminal velocity derived ear- 
lier (equation [41} . 

In a steady state, the solution of this equation implies 



Poq = pdo exp 



ie_ 

HI 



(69) 



where Hd is the "dust sphere" scale height, given by 

Hd = (2D4cd)'/' (70) 

Using equation |42] for sedimentation time tsed and c^ = 
GMfc/Ric, we obtain 



Rtc 



Sad Ef< 



A 



2 T,a X + a 



1/2 



(71) 
^ is the column depth 



where Efc = Mic/{-KRi^) ~ 10"^ g cm 
of gas in the first core f equation I14p. Recalling that Ea = 
paa, we see that for any grain sedimentation to take place 
(i.e., to have Hd <?C -Rfc), we need ad to be small and a to be 
large, e.g., macroscopic. For example, if ad = 0.001, a < A, 
and Pa ~ 1 g cm~^, grains larger than 1 cm are required. 



3.8 Destructive collisions 

Experiments show that the simple hit-and-stick picture for 
grain growth is modified when the relative velocity is larger 
than a few metres per second (jBlum fc Wurmll2008l ). Above 
these speeds grains can fragment or stick only partially. Fur- 
ther growth-reducing processes such as cratering or partial 
fragmentation can occur. To estimate the potential impor- 
tance of this, we recall that the terminal velocity of a grain 
of size a is 



47rGEa-R . . -1 R ^,,s-3/2 

Used = = 4.4 m s a -fr-T{t) ' 

Xs Hfc 



(72) 



where a <g A is in cm, and we assumed Mfc — 0.01 M0. Here 
the time-dependent dimensionless function T{t) (equation 
I30p accounts for the fact that the first cores contract with 
time, so that Rtc/cs oc pfc(i)"^/^ oc T(t)~^^^. 

From the above equation we see that growth of cm-sized 
grains may stall as their sedimentation velocity exceeds sev- 
eral metres per second. For a qualitative estimate of the 
effects of this, we can further assume that grain growth sat- 
urates at a size a such that Msed = Wmax, where Wmax is a 
few metres per second. In this case the grain sedimentation 
time scale is modified to 
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the addition of the radiative diffusion cooling: 
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Figure 4. Mean temperature (solid curves) of the first core of 
mass Mfc = lOMj as a function of time for a = 1 (red) and 
= 2 (black). The dotted curves show the virial temperature of 
the first cores in the analytical model presented in i|3.3l for same 
cases. Note that the analytical curves are within ~ 20 — 30% of 
the numerical calculation. 



tfii 



Ric/Vu 



3 X 10* yr 



_Rfc 3ms 
3 X IQia v^^^ 



(73) 



This "fragmentation-limited" sedimentation is slow, but not 
prohibitively slow compared with the grain growth time and 
the vaporisation times calculated earlier. The qualitative 
indication is hence that fragmentation of grains may slow 
down but not stop grain sedimentation inside the first cores. 



4 SET UP OF NUMERICAL MODELS 

Having learned about the problem perhaps as much as pos- 
sible within a simple analytical approach, we now turn to 
a spherically symmetric radiation hydrodynamics code to 
follow the evolution of the first cores and grain sedimenta- 
tion inside them. We use Lagrangian (gas) mass coordinates 
and the classical radiative diffusion approximation, which 
is appropriate given that first cores are very optically thick 
except for their atmospheres, which are of little interest to 
us here. 



4.1 No grain sedimentation case: gas only 
equations 

First of all we shall present evolution of the first core neglect- 
ing dust growth and sedimentation. Besides being interest- 
ing in itself, such a calculation is useful to contrast with the 
analytical solutions we obtained in i]3.3l The equations be- 
ing solved are the standard hydrodynamical equation with 



du 



_ = -A-kR 



de 



-47rP 



dP GMtotjR) 
dM 
d{R^u) 



dM 



i?2 

- 47r 



+ flg-d , 



d[R^F,, 
dET 
dR 



(74) 
(75) 
(76) 



where u is gas velocity, Mtot{R) ~ M{R) + Mg{R) is the 
total (gas -I- grains) mass enclosed within R, ag_d is the 
acceleration due to the gas-dust drag, to be discussed below, 
P is gas pressure, Frad is the classical radiation fiux, given 
by 



4a,adcr^ , p2 dT 

^-'^ = -^;<rr'"^ dM 



(77) 



where Orad is the radiation constant. The gas density is de- 
fined by p~^ = AnR^dR/dM. The numerical time integra- 
tion proc edure is based on the Lagra ngian scheme "Ihl" pre- 
sented in iBodenheimer et al.l (|2007l ). The scheme uses arti- 
ficial viscosity to capture shocks. We initialise the first cores 
as polytropic spheres of a given gas mass Afc as described 
in ^m 

Most of the simulations here do not include the convec- 
tive energy transfer since we did not expect it to be impor- 
tant: the optical depth of the first cores is not very high, 
so radiative cooling is quite efficient until the gas clump 
contracts significantly. This was confirmed in some of the 
simulations that did include the convective energy transfer. 
Furthermore, all of the simulations in paper II include the 
convective flux, and also the resulting convective mixing of 
the grains. We found that convection is never dominant be- 
fore the collapse of the central grain concentration into a 
massive solid core. Once that occurs, the release of gravita- 
tional energy from the core significantly changes the entropy 
profile near the solid core. This change drives very strong, 
in fact supersonic, convective motions in the region near the 
core. This strongly affects grain dynamics but after the solid 
core has already formed. 

The boundary conditions at _R = are it = 0, i\f (0) — 0, 
dP/dM = and Frad = 0. At the outer boundary we require 
continuity of gas density, velocity and the radiation flux. In 
this section we consider the case of dust tightly bound to the 
gas, in which case one can simply neglect the dust, setting 
ag_d to zero. 

Figure [4] shows the mean temperature in the cloud as a 
function of time for the flrst core of mass Mfc = O.OlAffc for 
the two opacity cases, a — 1 and a — 2. The dimensionless 
opacity coefficient fc* = 1 and the ambient gas temperature 
Ti — 1 for all the curves in the flgure. The slight jumps 
in the solid curves are due to initial oscillations of the first 
cores, which quickly decay away due to artiflcial viscosity 
(the numerical representation of the polytropic sphere initial 
conditions results in small perturbations). We also plot the 
analytical solutions (dotted lines) for the same clouds and 
opacity cases as obtained in ^3.31 Notably, the analytical 
and numerical solutions do not deviate from each other by 
more than a few tens of percent, which is more than can be 
expected from the order of magnitude analytical approach. 
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4.2 Inclusion of grain physics 

4-2.1 Grain dynamics 

We now turn our attention to the more complex case in 
which grain dynamics is taken into account. We use grains 
of the same size a for the whole cloud, where a is a free 
parameter. A more complex treatment is possible but is de- 
layed until a future paper. 

The grains are treated as a second fluid. We use the 
Lagrangian radial mesh provided by the gas coordinate to 
follow the evolution of dust. This is natural especially for 
small size grains that are strongly bound to the gas, as the 
relative gas-grain velocity is small, and the grains essentially 
move with the gas. The dust is described by a density in a 
given radial gas shell, and is allowed to slip through gas shells 
(see below). The grain initial density distribution follows 
that of gas, scaled down by the factor fg. 

The drag force acting on the grains/dust due to friction 
with the gas, (dua/dt) drag, is given by equation |39] Due to 
Newton's second law, the back reaction on the gas is 



'du\ 
IT I " ~^3 

\ aZ J drag 



dUa 

~dt 



(78) 



where Ua is grain velocity, pg is dust density (not to be con- 
fused with grain material density, pa\ and u is gas velocity. 
Accordingly, the acceleration ag_d in equation [75] is simply 

(du/dt)drag. 

The grains themselves are under the influence of gravity, 
drag force from the gas, and also undergo turbulent diffu- 
sion. From eauation l68l this can be accounted for by setting 
the grain velocity to 
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(79) 



where ui.^d is the grain sedimentation velocity, and the last 
term is the grain velocity generated by turbulent mixing. 
The terminal velocity approach that is used to derive equa- 
tion [75] is sufficiently accurate for small grains. For larger 
grains, the sedimentation velocity (equation I4ip may be 
larger than the local free fall velocity. To correct for this. 



Usad = mm[Wsed,"ff] 



(80) 



where mh = y/2GMtot{R)/r is the local free fall velocity. We 
tested this approach against the exact integration of grain 
radial time-dependent equation of motion, i.e., not assum- 
ing a terminal velocity. The differences are minor, and the 
terminal velocity approach is more numerically stable for 
smaller grains and does not require excessively small time 
steps. Therefore we pick the equation [TS] as the superior ap- 
proach for grain dynamics simulations here. 

Note that mass Mg{R), initially exactly equal to 
fgM{R), evolves with time separately from M{R), as grains 
are allowed to slip through the gas, and from one gas ra- 
dial mass shell into another. The mass continuity equation 
for grains inside a gas mass shell of index i, with inner and 
outer radii of Ri and Ri+i, respectively, is 



= - [R'^Psi-^a-u)],_^^+[R^Pg(Ua-u)], ,(81) 



1 dAMg 
4^ df 

where indeces i + 1 and i refer to the inner and outer 
boundary of the zone; e.g., pg^i is density of dust in zone 



i, AMg^i = (47r/3)pd,i(^?+i ~ Ri) is the grain mass inside 
zone i. 

There is a slight uncertainty in choosing the boundary 
conditions for grains at _R — >■ radius, that is in the very 
first gas mass zone. In contrast to the gaseous component, 
grains are not supported by pressure effects there, thus they 
can sediment and form a phase below even the very first 
mass zone of the gas. On the other hand, turbulence may 
suspend the grains in the fluid, curtailing the sedimentation 
into R — 0. Our goal here is to be conservative in our cal- 
culations of grain sedimentation. Therefore we assume that 
the turbulence keeps the grains suspended in the very flrst 
gas zone as long as the gas dominates the grains there by 
mass. Accordingly, the boundary condition for grain velocity 
is set to Ua{R = 0) = 0, as for the gas, for all the numeri- 
cal tests below. The simulations are stopped when the grain 
density in the inner zone reaches that of the gas density. 
One would expect that when the grain density exceeds that 
of the gas, the feedback that the grains impose on the tur- 
bulent motions of the gas becomes substantial, and hence 
the turbulent mixing support of grains should ease off. Pre- 
sumably a high densit y core comp osed of heavy elements is 
formed at this point (JBossJl 19981 ). A special treatment for 
the inner boundary condition is needed when this happens. 
We defer a study of the core formation process until a future 
publication. 



4-2.2 Grain growth and vaporisation 

We continue to consider only one size for the grains here, 
but in the interest of adding more realism to our models we 
allow that size to vary - grow by sticking with other grains, 
vaporise if gas enveloping the grain is too hot, and limit the 
grain growth if the differential grain velocity is too large. 

Our grain growth model follows the analytical prescrip- 
tion of H3.41 but with the addition of Brownian motions of 
small grains: 



da 



Pg 

4.pa 



{Uscd + Ubr) 



(82) 



where Used is the sedimentation velocity, and Ubr is the 
Brownian motion velocity of the smallest grains that 
dominate grain relative velo cities at small values of Usad 
(JDuUemond fc Dominikll200g ). The Brownian motion veloc- 
ity is a strong function of the smallest grains size, and we 
treat it as a free parameter of the model. The signs ( and ) 
signify grain-mass averaging of the quantities over the flrst 
core. This approach is a necessity given that the full radial 
and particle size distribution function treatment is beyond 
the scope of our initial study. 

Grain fragmentation or cratering may result from grain 
collisions with inter-grain velocities exceeding a character- 
istic value, which depend s on the size of the g rains partici- 
pating in the interactions (JBlum fc Wurml2008h . We explore 
these effects in a very simple fashion, introducing "maximum 
velocity" , iimax , above which grain growth is signiflcantly re- 
duced compared with equation 1821 In particular, we modify 
the above equation to 



da 



Pa 
4.pa 



{Uscd + Mbr) 



+ Wscd 



(83) 
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With this prescription, the da/dt — )■ as u^cd -^ oo. 

Finally, the grains are allowed to be vaporised if the gas 
temperature is lar g e. Th e vaporisation rate is taken from 
iHelled fc Schuberj (|2008l ). In practice this prescription im- 
plies vaporisation temperature between T — 1200 K and T 
= 1500 K, depending on grain size and relevant time scales. 
Vaporisation is not important for any models that do not 
reach 1200 K. 



5 NUMERICAL RESULTS 

5.1 Constant grain size, no turbulent mixing 

We start with the simplest, constant grain size case, setting 
a = 10 cm. We neglect the turbulent grain stirring, i.e., 
fixing ad = 0. The calculation assumes the grain material 
density of pa = 1 g cm^'^, initial mass fraction of fg — 0.005, 
opacity coefficient kq = 0.01, opacity power-law index a = 2, 
and the first core's mass of Mfc = lOMj. 

Figure [5] presents several snapshots showing the evolu- 
tion of the gas and the grain components. The snapshots 
correspond to times t ~ 360, 740, 1100 and 1500 years. The 
first inference from the Figure is that the constant density 
approximations for the first core itself and for the contract- 
ing grain sphere or "cluster" are reasonably good. The ana- 
lytical estimate of the settling time for this case is iaet = 500 
years (equation [42} . Further, equation 1561 predicts that the 
grains will collect into a central sphere the mass density 
of which equals that of the gas at time t — 880 years (for 
fg = 0.005). In the simulations, this occurs at time almost 
twice as long. This level of accuracy of the analytical esti- 
mates is nevertheless acceptable to us here. The simulations 
show a somewhat more complex evolution than the simple 
"top hat" profile. 

The upper right panel of the figure shows the radiation 
fiux as a function of radius within the first core. Note that 
in the last snapshot the heat flux increases significantly in 
the inner ~ 10% of the first core. This is a signature of the 
adiabatic contraction heating imposed by the ever increas- 
ing grain density in that region. While until this point the 
effects of grain sedimentation were hardly felt by the gaseous 
component anywhere inside the core, the "inverse drag" - 
grain drag on gas is now significant in the inner region. The 
gas is however stable against further collapse there as can 
be seen from velocity curves (lower right panel): the gas and 
the grain phases continue to separate out, e.g., move with 
different velocities. The reason for which the gas phase re- 
mains gravitationally stable is that it takes a relatively small 
amount of heating (the bump in the flux curve in the upper 
right panel) to set up a sufficient pressure gradient to oppose 
the collapse. 

The last snapshot of the figure corresponds to a stage 
between the grain-dominated one f ^3.6.2|) and the "bound 
grain cluster" ( H3.6.3p . Soon after the last snapshot shown 
in the figure the grain distribution undergoes a dynamical 
collapse which we shall study in a future paper. 
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Figure 6. Density (solid for gas, red dotted for grains) and tem- 
perature distributions at three different times for the same calcu- 
lation as in Figure [5] but now with the turbulent mixing coeffi- 
cient a^ = 10~^. The snapshots are for times t = 0, 1100 and 2200 
years. Note that the grain distribution concentrates towards the 
centre by a small amount only before coming to an equilibrium 
in which sedimentation is balanced by turbulent mixing. 



5.2 Constant grain size, turbulent mixing 

5.2.1 Strong turbulence 

We now add turbulent mixing, still keeping the size of the 
grains constant at a = 10 cm. Figure [6] shows the density 
and temperature distributions for the same case computed 
in i]5.1l but for the turbulent diffusion coefficient Od = 0.01. 
It is obvious from the snapshot sequence that the dust dis- 
tribution quickly adjusts to an equilibrium shape in which 
the rate of gravitational settling is offset by turbulent dif- 
fusion. The contraction of the density distribution between 
the last two snapshots is entirely due to the contraction of 
the gas distribution, which occurs on the radiative cooling 
time scale of about 10^ years. Note that no dense grain clus- 
ter forms in this case no matter how long we were to follow 
the calculation. After the first core has contracted enough to 
reach T ^ 2000 K, it would undergo the well know n second 
collapse (|Larsonlll969l : lMasunaga fc Inutsukall200d ). Appar- 
ently, the effects of grain sedimentation for these particular 
parameters are completely negligible. 
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Figure 5. Grain sedimentation for a first core of mass Affj. = lOAIj and a fixed grain size of a = 10 cm. Tlie panels show snapsiiots 
of gas density, temperature, radiation fiux and velocity, as labelled, in solid curves. The dotted and dashed red curves show the grain 
density and velocity, respectively. The snapshots are for times between t = 360 and 1500 years. 



5.2.2 Weaker turbulence 

Figure [7] shows a calculation identical to that presented in 
tj5.2.1l and in Figure [S] but now for a weaker turbulent mix- 
ing, with Qd = 10""^. Once again an equilibrium is quickly 
set up, but a much more concentrated one. It is interest- 
ing to note that this equilibrium implies a strong metalicity 
gradient within the gas cloud, with the most metal rich gas 
being in the central regions. 



5.2.3 Weaker turbulence and larger grains 

Our final experiment with fixed size grains is to repeat the 
calculation of i|5.2.2l but with a larger size of grains, a = 30 
cm. Figure [5] demonstrates a very rapid sedimentation of 
grains in this case, with gravitational collapse of the "grain 
cluster" f ^3.6.3p occurring by about 500 years. As expected, 
larger grains are affected significantly less by the turbulent 
mixing. 

The rapid grain sedimentation in this calculation yields 
a mean grain sedimentation velocity of about 15 m s^^ in 
this model. Such a rapid sedimentation would undoubtedly 
lead to shattering of some of the grains, so that this calcu- 
lation is somewhat unrealistic, but it goes to show that if 



grains can grow to decimetre sizes, then they can probably 
overcome the effects of turbulent mixing. 



5.3 Full models 

5.3.1 Simulations Table 

We shall now explore our "full" model that includes the 
grain growth prescriptions ( i|4.2.2[) . Selected runs that we 
performed and discuss below are listed in Table [T] Some of 
the parameters of the models have very little bearing on the 
results, except near the boundaries separating qualitatively 
different regimes. These parameters are not listed in the Ta- 
ble. In particular, the results are insensitive to the initial 
size of the grain, ao, and the Brownian motion velocity, Ubr, 
as long as the latter is large enough (most of the runs were 
performed with Ubr = 5 cm s~^). This is because, start- 
ing from /im-size grains, one finds that grains quickly reach 
sizes of hundreds of nm anyway (see also more detailed cal- 
culations, showing a similar initial rapid grain growth, by 
IDuUemond fc Domin ik 2005). Thus ao is set to lOOfim for 
most tests below. 

Given the large differences in the cooling times ( H3.ip . 
the form of the opacity law predictably turns out to be 
quite important. The names for the runs performed with the 
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Figure 7. Same as Figure |6] but for a smaller turbulent mixing 
coefficient, a^j = 10~^. The dust is concentrated towards the 
centre much more than in the previous calculation, but still hovers 
at an equilibrium state. 



power-la-w index a = 1 starts with a letter "S" for "small 
opacity", -whereas names for a = 2 runs start with "L" for 
"large opacity". These names are not to be taken literally 
as there is also the opacity coefficient, ko, in front of the 
opacity law we use here (equation [T]) . kq is also varied in 
some models. 

The last three columns in table [l] list several output 
quantities of the simulations, igrain , Tc and pc are the time 
at which the self-gravitating phase of the grain sphere evo- 
lution is reached (if ever), and the central temperature and 
the density of gas, respectively, at that time. Grains were 
vaporised before they could make a massive self-gravitating 
sphere in the centre in the runs in Table 1 for which no 
values of these three parameters are given. 

5.3.2 Turbulent mixing versus grain growth 

We first discuss runs LI, L2 and L3, the three models that 
are identical to each other except for the value of the tur- 
bulent mixing coefficient, which was set to aa ~ 10"'*, 10"'^ 
and 10~^, respectively. There is very little difference in the 
evolution of the models. This may appear surprising given 
the all-important role the turbulence played in the simula- 
tions with the fixed grain size. The "paradox" is resolved 
by realising that the runs always start with the turbulence- 



Figure 8. Same as Figure [7] but for larger grains, a = 30 cm. 
The times of snapshots are t = 0, 220, and 440 years. Soon there- 
after, the inner grain-dominated part of the first core undergoes 
a gravitational collapse similar to the a^j = case studied earlier 
(Figures [Sjl . 



dominated regime (small grains) and end up in the negligi- 
ble turbulence one (large grains) , and only the time and the 
boundary between these two vary in L1-L3. 

Figure |9] shows the grain size, a, versus time for these 
three models. The curves are terminated at the point when 
the central radial bin is dominated by the grains. Initially the 
curves in the figure are nearly identical to each other. This 
can be understood by realising that at a small a, sedimenta- 
tion speeds are low, so that turbulent mixing easily balances 
sedimentation. A quasi-steady state results in which grain 
sedimentation is neutralised by turbulent mixing. 

Figure [9] shows that grain growth is slightly faster for 
higher levels of turbulence. We explain it by the fact that 
the grain distribution is slightly more extended in this case, 
hence the sedimentation velocity, proportional to R (see 
equation I4ip . is larger, allowing the grains to sweep up 
smaller ones at a faster rate and hence grow faster. 

The process of grain growth accelerates with time. At 
early times this is due to the fact that Used oc a, thus 
da/dt ex. a (see equation I82|l . At later times, grains become 
over-abundant in the inner regions, i.e., pg/p > fg there. 
This speeds up the rate at which the larger grains can sweep 
up the smaller grains even further. This non-linear stage of 
grain growth is naturally reached faster by the lower turbu- 
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Figure 9. Grain size versus time for models described in i|5.3l 
The solid black, the red dotted and the blue dashed curves are 
for a^i = 10~*, 10"'^ and 10"'^, respectively. 



lence models. A corollary of this is that the higher turbu- 
lence runs reach the grain self-gravitating phase with higher 
temperature and density. 

Nevertheless, it is notable how qualitatively similar the 
three simulations actually are. One difference is the time 
when the central region becomes grain-dominated. It is i = 
7.6, 8.4 and 11 thousand years for a^ = 10"*, 10"^ and 10"^, 
respectively. The later the grain-dominated phase reached, 
the hotter is the central region of the first core, and the 
denser it is (see Table 1). 

The main conclusion we can draw from these experi- 
ments is that the turbulent diffusion is only likely to de- 
lay rather than stop grain sedimentation if grain growth 
can proceed via the hit-and-stick process faster than the gas 
heats up to vaporise the grains. 



5.3.3 Delayed sedimentation due to grain fragmentation 

To explore sensitivity of our grain growth model to grain 
fragmentation due to high speed collisions, we ran several 
models where the maximum velocity parameter Umax was 
reduced from the nominal value of 10 m s~^ to 1 m s~^ 
(see Table [T]). In particular, simulation L2a is identical to 
L2 except for the smaller value of Umax. Qualitatively the 
two runs are similar, but the grain-dominated cluster phase 
is reached at time almost twice as long in L2a than in L2. 
This is caused by the subdued grain growth rate in the run 
L2a compared with L2. 

Concluding, as envisaged earlier in H3.8I grain sedimen- 
tation may still proceed reasonably quickly at realistically 
small ('^m s~^) values of Umax- 



5.4 Vaporisation of grains in rapidly cooling first 
cores 

In the context of our model of isolated first cores, the only 
robust way to stop grain sedimentation is to prevent their 
growth. This occurs when the first cores become too hot, so 
that grains vaporise. We find this behaviour if (a) the cores 
are too hot (too massive) to begin with; (b) the opacity 
is low so that the cores cool very quickly (thus leading to 
hotter internal temperatures); (c) the grain content of the 
core is low, yielding too long grain growth times. 

Figure [10] shows the grain size and temperature as a 
function of time for three runs that yielded no grain sedi- 
mentation. The solid black curve is the simulation S3 (cf. 
Table 1) - a low opacity, kq — 0.001, and low grain mass 
fraction case, fg = 1.2 x 10~^, with the first core mass of 
Mfc = 5Mj. The opacity law exponent is a = 1. Both con- 
ditions (b) and (c) are the case for this particular run. The 
grains manage to reach the size of about 1 cm before they 
are vaporised. The central gas temperature is almost 1400 
K at that time. Grain vaporisation actually starts earlier, at 
Tc ~ 1300 K or so, but does not immediately prevail over 
the process of grain growth. 

The other two runs shown in figure [TOl have much more 
massive cores, Mfc — 25 Mj (S7 in Table 1; dotted red curve) 
and Mfc = 20Mj (S6; dashed blue). The opacity coefficient 
fco = 0.05 and 0.01, for the runs, respectively, and the grain 
mass fraction is fg = 0.01 for both. The cores cool quickly 
nevertheless as the opacity power law index a = 1, hence 
the short cooling times (equation 1211) . The most massive 
first core considered here (blue curve) is almost too hot to 
allow grain growth to begin with (case "a" above). As the 
cooling time is short, the core undergoes a quick and nearly 
dynamical evolution in the beginning, driving hydrodynami- 
cal oscillations that allow some early grain growth. However 
the latter is quickly turned into grain evaporation as the 
core settles into its cooling evolution. 



5.5 Importance of grain fraction (metalicity) 

Amongst the parameters of the simulations that we varied 
is the first core's initial grain mass fraction, fg (Table 1). 
A higher grain content obviously favors grain growth. For 
example, the runs L4 and L5 are identical except for the 
grain mass fraction, fg = 0.005 and 0.02, respectively. Not 
surprisingly, the grain-dominated phase is reached far faster 
in L5 than in L2. Therefore, we can expect that metalicity of 
the first core is one of a key parameters determining whether 
grain sedimentation takes place or not. 



6 DISCUSSION 

6.1 Summary of the results 

In this paper we explored the possibility that grains may 
sediment inside isolated or slowly accreting first cores, which 
are the gaseous condensations with mass between about 5 
and 50 Jupiter masses. While many parameters (such as 
the turbulent viscosity parameter) determine the end result, 
the main requirement for sedimentation to take place is a 
rapid enough grain growth. Based on numerical experiments, 
sedimentation occurs as long as grains can grow to a few 
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Table 1. Parameters of the simulations. The mass of the first core, Mf^, is in Jupiter's masses; fco is the opacity coefficient in the opacity 
law given by equation [71 o is the power law index in that law; fg is the grain mass fraction of the core; a^ is the turbulent mixing 
coefficient; fmax [m s~^], the maximum velocity before grain-grain collisions stifle grain growth (equation I82I I; tg,ain is the time at which 
the innermost gaseous cell becomes dominated by grains, by mass, in units of 10^ yrs; and T^ and pc are the gas temperature and density 
there at the time, in cgs units. 



ID 


Mfc 


fco 


a 


h 


Qd 


^max 


*grain (103 yrs) 


Tc 


Pc 


SI 


fO 


0.01 




0.005 


5 X 10-2 


10 


2.7 


1400 


1.4 X 10-* 


S2 


fO 


0.001 




0.005 


5 X 10-2 


- 


- 


- 




S3 


5 


0.001 




0.0012 


10-3 


1 


- 


- 


- 


S4 


5 


0.003 




0.0012 


10-3 


1 


3.3 


1270 


1.5 X 10-* 


S6 


20 


0.01 




0.02 


10-3 


1 


- 


- 


- 


S7 


25 


0.05 




0.02 


10-3 


1 


- 


- 


- 


Lf 


fO 


0.01 


2 


0.005 


10-" 


10 


7.8 


695 


7.5 X 10-1" 


L2 


fO 


0.01 


2 


0.005 


10-3 


10 


8.4 


710 


8 X 10-10 


L3 


fO 


0.01 


2 


0.005 


10-2 


10 


11.0 


750 


10-9 


L2a 


fO 


0.01 


2 


0.005 


10-3 


1 


13.4 


790 


1.2 X 10-3 


L4 


fO 


0.001 


2 


0.005 


5 X 10-2 


10 


1.9 


1190 


4 X lO-'' 


L5 


fO 


0.001 


2 


0.02 


5 X 10-2 


10 


0.69 


1130 


3 X lO-'' 


L6 


fO 


0.1 


2 


0.005 


10-3 


10 


16.0 


360 


1.3 X 10-"' 


L7 


40 


0.1 


2 


0.005 


10-3 


10 


15.3 


1100 


2.1 X 10-1" 


L8 


55 


0.1 


2 


0.005 


10-3 


10 


- 


- 


- 



cm to a few tens of cm before vaporisation temperature of 
T ~ 1400 K is reached. This is expressed mathematically as 



tgr *^ ''V 



(84) 



where igr is the grain growth time scale (equation I49|l . and 
tvap is the vaporisation time of the grains (eauation l33p . Tur- 
bulent mixing and grain shattering by collisions may slow 
down grain growth compared with the analytical model (in 
which case tgr must be obtained numerically). 

While the outcome depends on several parameters - 
mass of the first core, Mfc, ambient gas temperature, Tlnit, 
opacity law in the first core, the initial grain mass fraction, 
fg, turbulent mixing, etc., we shall here give only one ex- 
ample of the implications of the equation 1841 In particular, 
we shall fix fg — 0.005, Od = 10-^, Umax = 10 m s"i, and 
Tinit = 10 K. 

With these parameters fixed, we can ask the following 
question: At a given opacity coefficient ko, what is the range 
of the first core masses that could support the dust sedimen- 
tation? 

Figure [11] answers this question for the two opacity 
power-law indexes that we studied here, a — 1 (black dotted 
curve) and a = 2 (red dashed one). The lower (solid) line 
gives the minimum mass of the first core which we obtained 
by requiring that the mean first core density (equation llip is 
at least equal that at which the collapsing gas switches to the 
adiabatic behaviour (equation [8]). This curve is independent 
of a for Ti = 1. The curves dividing the parameter space 
are obtained by solving the equation [84] with our analyti- 
cal model for grain growth and the first core contraction. 



The asterisks are the maximum first core mass for which 
grain sedimentation occurs found numerically for the a — 2 
case. They are encouragingly close to the analytical curve. 
The fact that they are slightly below is probably explained 
by the extra delay in the grain growth that the turbulent 
mixing and the maximum velocity «max impose in numer- 
ical simulations compared with the analytical models that 
do not account for these processes. 

This limited survey of parameter space demonstrates 
that there range in Affc in which grain sedimentation is pos- 
sible is about a factor of 10 wide for a = 2, and shifts to 
higher masses for higher opacities. This is to be expected, 
as higher opacity implies longer cooling times. 

The Q = 1 case (black dotted line) shows a narrower 
window of opportunity for dust sedimentation, which is 
again a cooling effect. As these cores cool quicker, there is 
less time for dust growth. 

6.2 Comparison to previous work 

Our analytical mode l and numer ical simulations confirm the 
suggestion made by iBosd (|l998l ) that grains can grow and 
sediment inside giant embryo, although started from a much 
less dense initial configuration. While a fuller comparison 
is to be made in paper II where we continue the calcula- 
tions to the point of the solid core formation, we can al- 
ready see that this lower density and temperature initial 
configurati on is key to address the c r iticism levelled on the 
JBossI (Il998l) model bvl\yuchterl et all (120001), late r confirmed 
bv lHelled etHI (|2008i ): lHelled fc SchuberT(|2008r ). These au- 
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Figure 10. Grain size versus time (upper panel) and the central 
gas temperature (lower panel) for three cases that led to grain va- 
porisation rather than growth: S3 (solid black), S6 (dashed blue) 
and S7 (red dotted). In all three cases the first core cools radia- 
tively, contracts, and becomes too hot to allow grain growth. 



thors found that their giant planet models were convective, 
and that convective grain mixing has significantly slowed 
down the core growth. 

Since the initial column depth of the embryos is much 
lower in our models, the radiative cooling remains strongly 
dominant over convective cooling for ~ 10 — IC years, 
the time it takes (see figure 2) for the embryos to con- 
tract to densities and temperatures considered by the au- 
thors mentioned ab ove. These findings echo the results by 
iBodenheimen (| 19741 ) who studied early contraction of a 1 
Jupiter mass gas cloud and concluded that the planet re- 
mained radiative during the first ~ 10'^ years. While his 
model does not include dust sedimentation, it accounts for 
H2 molecules dissasociation, hydrogen ionisation, and much 
more (that we do not include here). His calculations also 
used similarly cool initial conditions (43 K was the initial 
temperature of his model planet). 

Therefor e it appears that most of the differences be - 
tween us and lHelled et all (|2008l ): lHelled fc Schubert (|2008l ) 
can be traced to the significantly different initial condi- 
tions. In terms of applications to plane t formation fi e ld, th e 
differences in our r esults and that of iHelled et al.l (|200a ): 
iHelled fc Schubert] (|2008l ) are striking. These authors found 
that a 5 Jupiter mass embryo cannot yield any dust sed- 
imentation as it is too hot to begin with, e.g., the initial 
temperature above the grain vaporisation temperature. In 
contrast, we find that a 5 Jupiter mass embryo provides an 



excellent environment for dust growth: the gas temperature 
is initially less than 100 K and stays below 1400 K for as 
long as 10* to more than 10^ years, depending on the opacity 
(cf. figure 1). 

In paper II we do find that convection becomes very 
important in the inner region of the embryo once a mas- 
sive solid core forms, as this releases a significant amount 
of heat. It is also interesting to note that while our initial 
condition is significantly different, the inferred grain growth 
and sedimentation time scales are from 10'^ to 10* y ears, in 
a broad agreement with the estimate bv lBosg (| 19981 ). 

On the other h and, in a qualitati v e agr eement with 
IHelled et al] (|2008l ): IHelled fc SchubertI (|2008l ). the condi- 
tions for dust growth found to deteriorate with an increas- 
ing mass of the gas clumps. Depending on opacity, grains 
cannot sediment for giant embryo masses higher than ~ 10 
to a few tens Jupiter masses. 



6.3 Astrophysical implications 

Our calculations h ere and in p aper II lend support to the 
ideas presented by iBosd (|l99g ). except that grain sedimen- 
tation process must start earlier in the life of the proto- 
planet, while it can still be considered the first core. Lower 
initial temperatures and rather long cooling times of the gi- 
ant embryo in the latter stage are key for a successful grain 
sedi mentation outcom e . 

IBoss et al.l (|2002l ) argued that giant embryos may yield 
not only giant planets but also giant icy planets with cores 
if the metal poor envelopes of the embryos are removed. 
Irradiation by a nearby star was suggested to accomplish 
this. While this is certainly possible, it may be not very 
probable as OB stars are rare. 

The present paper is the base for paper III, where we ar- 
gue that giant embryos may migrate inward due to the gravi- 
tational torques from the disc. Estimates show that embryos 
may migrate to the distance of several AU to the parent star 
in several x 10* years, typically. As cooling of embryos slows 
down as they age, the embryos must be eventually disrupted 
by the tidal forces or by the heating due to irradiation from 
the parent star. This opens up an exciting possibility that 
all planets may be formed by such a modified version of 
gravitational instability model. 



7 CONCLUSIONS 

In this paper we considered the grain growth and sedimen- 
tation process inside the giant planet embryos, which we ar- 
gued must be nearly identical in properties to the first cores. 
These are the first gaseous hydrostatic condensations from 
which stars may form. We found that grains can indeed grow 
and sediment to the centre of the gas cloud provided that 
the gas remains cool enough (temperature below ~ 1400 K) 
for the time it takes the grains to reach a few cm size. The 
efficiency of the dust sedimentation process and the final 
mass of the solid core are strong functions of opacity and 
other parameters of the problein. 

We suggested that astrophysical applications of these 
results may be in the field of planet formation, where giant 
planet embryos may serve as birth places for all types of 
planets if these embryos migrate inward and get tidally or 
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Figure 11. Minimum first core mass (solid curve) and maximum 
first core mass for grain sedimentation versus the opacity coef- 
ficient Kg for o = 1 (blaclt dotted line) and a = 2 (red dasfied 
line). No grain sedimentation is expected above these lines as 
grains get vaporised faster than they can sediment. See text for 
further detail. 
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